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Abstract : For interconnection of two or more areas in power 
system , frequency should be maintained within the scheduled 
value, which can be achieved by employing one of the most 
prominent techniques called as Automatic Load Frequency 
Control (ALFC). In ALFC, the frequency can be controlled in 
three ways, namely Flat frequency regulation, Parallel 
frequency regulation and Flat tie-line loading . Among these 
controls, Parallel frequency regulation is commonly used 
method, because constant frequency can be maintained by 
equalizing the power generation with the power demand. The 
main aim is to reduce the oscillations or damping in the 
frequency and also the changes in tie-line power. In this paper, 
PI and PID controllers are used to improve both frequency and 
tie-line power responses for multi-area interconnected power 
system (which consists of Hydro-thermal generating system ) 
and is developed based on the model predictive control. This 
paper presents a solution procedure to Load Frequency Control 
(LFC) of interconnected Hydro-Thermal system with a novel 
controller (PSO (Particle Swarm Optimization)). 

Keywords: Automatic Load Frequency Control, Conventional 
Controllers, PSO technique, Area Control Error (ACE). 

I. Introduction 

Generally, power system consists of three parameters which 
shall be within the limits for successful operation i.e. Frequency, 
Voltage and Load angle, among these frequency parameter plays 
vital role [2] . Many different power frequencies were used in the 
19th century. Very early isolated ac generating schemes used 
arbitrary frequencies based on convenience for steam 
engine, water turbine and electrical generator design. Frequencies 
between \6 2 A Hz and 133 1 /3 Hz were used on different systems. 
The main purpose of a power engineer is to provide power to the 
consumers reliably and economically with a better quality [4-5]. 
The frequency and tie-line power should be kept within the limits 
by equalizing the power generation at the generating end and the 
power consumption at the load end, because there are two points 
available throughout the power system for keeping the frequency 
within the limits, one is at the generating end and the other is at 
the load end. Load frequency control problem arises when 
individual generation areas are interconnected by transmission 
lines called as tie-lines. 

The large-scale power systems are liable to performance 
deterioration due to the presence of sudden small load 
perturbations [1], parameter uncertainties, structural variations, 
etc. Due to this, modern control aspects are extremely important 



in Load Frequency Control (LFC) design of power systems [3]. 
The main objectives of Load frequency control are i) to 



Area-1 Area -2 




maintain constant frequency throughout the system, ii) to 
preserve the tie-line power at a scheduled level irrespective of 
load changes in any area, iii) to diminish Area Control Error 
(ACE) and iv) to minimize Peak time, Peak overshoot and 
Settling time. When interconnecting two or more areas/regions 
for providing continuous supply, the complexity of the system 
increases. Due to the increased complexity of modern power 
systems, advanced control methods were proposed in ALFC, e.g., 
Fuzzy logic controller, evolutionary based techniques, 
decentralized optimal control, pole placement technique, APSO 
(Accelerated Particle Swarm Optimization) and Quasi 
Decentralized Functional Observer (QDFOs) [14]. The most 
important features of two utility interconnected multi-area 
thermal system is that a) to buy or sell power with neighboring 
systems, b) to meet sudden requirement of power demand, c) to 
reduce the installed capacity, d) to possess higher reliability with 
better quality of supply. The Modeling and the simulation of two 
area Hydro-thermal interconnected power systems [6] are 
realized in MATLAB Software. 

This paper is organised as follows, section-I explains about the 
introduction of Load frequency control and how the situation 
exist in early years of electricity generation. Section-II briefly 
describes about the system. Section-Ill illustrates about different 
control strategies such as conventional controllers, PSO 
technique. Section-IV describes about mathematical modeling of 
each and every part of single area load frequency controller with I, 
PI and PID controller. Section -V gives you about two area 
hydro-thermal interconnected system using the specified 
controllers in section-III, and Section- VI shows you the results 
with comparison of important specifications such as Peak 
overshoot, settling time, steady state error, etc. of frequency of 
each area and tie-line power. 
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II. System Investigated 

The LFC of interconnected system consists of two generating 
areas of equal size; area-1 includes electrical governor hydro 
system and area-2 make up a reheat thermal system. Fig. 3 shows 
the FFC model of single area thermal system with integral 
controller and fig.4 shows the AFFC model with electrical 
governor in hydro area and single stage reheat turbine in thermal 

area. A frequency bias setting (3 j is considered in both hydro and 
thermal areas. MATFAB/SIMUFINK has been used, to obtain 
dynamic responses for Afj , Af 2 and AP tie _ line for 1% step load 

perturbation in either area. In the literature, the performance of 
the electrical governor was found quite superior than the 
mechanical governor [6] ; so that in this paper i included electrical 
governor only in hydro area. The optimal values of derivative, 
proportional and integral gains for the electrical governor, speed 
regulation droops i.e. R b R 2 , synchronizing coefficient T 12 , time 
constant and gains of hydro, thermal turbine and power system 
have been taken from the work of paper [15]. 

III. System Optimization using PSO 

A proportional-integral-derivative controller is a control 
loop feedback mechanism widely used in industrial control 
systems. It calculates an error value as the difference between a 
measured process variable and a desired set point. The response 
of the controller can be described in terms of the responsiveness 
of the controller to an error, the degree to which the 
controller overshoots the set point, and the degree of system 
oscillation. 

PSO is an optimization technique which provides an 
evolutionary based search. This search algorithm was introduced 
by Dr. Russ Eberhart and Dr. James Kennedy in 1995. It was 
inspired by the social behavior of bird flocking and fish schooling. 
The basic concept of this technique that the change in the 
velocity of each particle towards its P best and g best positions at 
each step. This means that each particle is try to modify its 
current position and velocity according to the distance between 
its current position & P best and the distance between its current 
position & g best . 

The parameters of controller in both the areas have been 
optimized using integral square error (ISE) technique. The 
objective function for this technique, to minimize frequency 
deviations in area-1 & 2 and tie-line power deviations, used is 

J = j” [(Af, ) 2 + (Af 2 ) 2 + (AP tie _ lme ) 2 ]dt (1) 

0 

Subjected to, 

■^il(min) ^ ^il ^ ^il(max) 

Where, Af x , Af 2 and AP tie _ _ Une are the incremental changes in 

frequency of area-1, area-2 and the incremental changes in tie- 
line power respectively. 

PSO allows us to intelligently search the solution space for a 
global best solution. In contrast to an exhaustive search, PSO 
significantly reduces the time to optimize response. It was 
originally developed by biologists to be able to mathematically 
describe the swarming of birds. 



The basic PSO algorithm is composed of two 
mathematical steps namely i) intelligently update the velocity of 
the particle and ii) update the position of the particles based on 
their velocity. The velocity of each particle is update according to 
the following equation, 

v- +l = wv k + Cri k ( Pl k - x k ) + sr k (p k - x k ) (2) 

Where, xf is the position of i th particle at k th iteration 
vf is the velocity of i th particle at k th iteration 

p A is the best solution observed by i th particle 
p is the best solution observed by any particle 

W is the inertia coefficient 
C is the cognitive coefficient 
S is the social coefficient 
r is the randomness factor 

In the above equation, the first term refers to inertia 
coefficient which controls how quickly a particle will change 
direction, the second one represents cognitive coefficient which 
controls the tendency of a particle to move towards the best 
solution observed by that particle and the last term expresses 
about social coefficient which controls the tendency of a particle 
to move towards the best solution observed by any of the 
particles. 




Fig.2 Flow chart of PSO Technique 
The position of each particle can be update using the following 
equation 



X k+1 =x k 



■ XV, 



k+1 



(3) 
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Where % is called as constriction term which is used to tighten 
up the group towards the end of optimization. 



IV. Model of Single Area Load Frequency Control 



The system modeling will be developed briefly by 
considering two assumptions. The first is that the incremental 
changes in power demand is neglected, the second is that control 
of real power and frequency, and control of reactive power and 
voltage, are decoupled and can be considered separately. 

A. Speed Governing System: 

AX e (S) = f { Af (S), AP C (S)} (4) 



Calculation of AX E (S) : 



AX e (S) 



Kq 

1 + ST g 



AP C (S)- 



AF(S) 

R 



Where, 

Static gain of the governor, 
Time constant of governor, 



K 0 = 

T g = 



K 2 K 3 

~k7 

i 

K 4 K 5 



and 



(5) 



Regulation of speed governor or droop value, r = A 

K, 



B. Turbine Model: 

There are two types steam turbines namely reheat turbine and 
without reheat turbine. 

Rating of high pressure turbine stage is Ot p.u. say 

Both the ratings of intermediate and low pressure turbine stages 

are 1 — Ot p.u. 

AP V (S) is Turbine input 

AP t (S) is Turbine output 
Total transfer function of turbine 



AP t (S) = AP T HP (S) + AP t ip LP 
K t 1 + cxST R 
1 T- ST t 1 + ST R 



(S) 

AP V (S) 



Where T R is reheating time constant, 

K t is static gain of turbine, 

T t is turbine time constant 

C. Generator model: 

Net accelerating torque, 

Let or T a = T mech - T elec 
= Ixa 

Angular momentum of machine, M = cox I 
Where, T mec h is mechanical input torque in N-m 
T elec is electrical output torque in N-m 
I is inertia constant or moment of inertia of machine 
Ot is angular acceleration in rad/sec 2 
CO is angular speed in rad/secs 
AP mech (S) - AP elec (S) = MSAco(S) 

D. Model of load: 

Change in frequency dependent load, AP Lf = D x Aco 

D = AP^ MW/Hz or (p.u. MW)/(p.u. Hz) 

Af 



( 6 ) 



(7) 

( 8 ) 



(9) 



Aco(S) = [AP mech (S) - AP L (S) - DAco(S)] 

Where D is the damping constant 
E. Combined model of Generator and load: 



G P (S) 



1 K p 

MS + D 1 + STp 



Where K P is power system gain = 1/D 

T P is power system time constant in secs =M/D 

AP C = -J KjAf (t)dt AP C (S) = - A-AF(S) 

s 



( 10 ) 



(ID 



( 12 ) 




Integral controller Governor Non-reheat turbine Power system 



Fig. 3 MATLAB/SIMULINK model of Single area thermal 
system with integral controller 

There are two loops in the ALFC namely primary and 
secondary loops. Primary loop comprises model of Fly ball speed 
governor, Non-reheat turbine and combined model of generator 
and load called as power system. Power system consists of both 
frequency dependent and independent loads, in this paper 
frequency dependent load only included. After modeling the 
primary ALFC loop there is always some error present in the 
output at the steady state. 

This error could be eliminated by introducing proportional 
and integral controller then the overall model can be called as 
secondary ALFC loop. This error can be reduced with ease using 
conventional controllers. For single area ALFC, with PID 
controller the response of frequency could be settled in less time 
when compared to I and PI controllers. 

V. Two Area Automatic Load Frequency Control 

The frequency should be constant whenever you interconnect 
two or more areas in AC power systems. For operating these 
areas simultaneously the variations in the tie-line power should 
be kept within limits. Care should be taken while designing Tie- 
line, since it has the capability of carrying over power demand by 
the neighboring areas. The modeling equations of tie-line power 
as follows: 

Model of Tie-line: 

Power flow from areal to area2 without disturbance 

P 12 = sin(8° — 8° ) (13) 

^12 

Net power flow from areal to area2 due to disturbance 

P i 2 = AA[sin( 5 ? -S°)]+ AA[cos( 5 ? - 6 °)(A 6 , - A 5 2 )] 

^12 ^12 

AP 12 = T 12 (A5 1 - A5 2 ) (14) 
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Where, 

T 12 is known as synchronizing coefficient x 12 = VlV2 -8°) 

X 12 



(15) 

By taking Laplace transform, 

ap 12 (S) = ^^[af 1 (S)-af 2 (S)] (16) 

The MATLAB/SIMULINK model of two area interconnected 
Hydro- thermal power system is shown in Fig. 6. In this the input 
to the controller is the sum of tie-line power and the product 
frequency bias parameter p and change in frequency Af. i.e. 

ACE! = AP 12 + P x Af x (17) 

ACE 2 = — AP 12 + P 2 Af 2 (18) 

Where and P 2 are frequency bias constants of area-1 & 2 



respectively and AP 12 is the change in tie-line power between the 



two areas. 

The speed of generator and thus frequency of the overall 
system with respect to changes in the demand could be achieved 
by choosing optimum values of frequency bias, regulation droop 
and the gains of particular controller. In this way the active power 
generation can be matched with the power demand consumed by 
the load. 




seconds, because in thermal and hydro power systems consists 
mechanical equipments which has more time constants when 
compared with the electrical devices. 

The conventional integral controller gains for the hydro area 
and thermal area are found to be K n =-1.2 and K I2 =-0.5 
respectively. Using conventional PI controller, optimum gains 
Kpi=-0.01, K n =-0.6 for hydro area and K P2 =-0.1, K I2 =-0.5 for 
thermal area. The optimum value of integral controller gains for 
the hydro area and thermal area are K n =-0.84 and K I2 =-0.98 
respectively, using ISE technique. 




Fig. 5 Comparison of I, PI and PID controller response of 
single area non-reheat thermal area 

The steady state error and time response specifications such as 
peak overshoot, peak time and settling time are improved and are 
clearly shown in Fig. 5. The comparison among these parameters 
in quantitatively has shown in table. 1 with I, PI and PID 
controllers. By investigating these controllers the peak overshoot 
is reduces from 0.52 to 0.09. Here the Settling time (the time 
taken by the response to reach steady state value), was considered 
for 2% tolerance band, is less for integral controller then for PI 
controller and more for PID controller. The variation in the delay, 
rise and peak times of these controllers are quite similar. 



Fig.4 MATLAB/SIMULINK model of two area hydro-thermal 
interconnected system 

As the number of interconnecting areas increases the 
complexity of the system will increases, so it is difficult to control 
the system with conventional controllers for the corresponding 
changes in the load demand. The power system researchers are 
going to adopt more and more optimization techniques to control 
the valve input to the turbine as fast as possible. 

VI. Result and Analysis 

The results of single and two area interconnected hydro - 
thermal power systems are clearly analyzed in the following: 
Here the time taken to reach steady state value is in the order of 



The responses of Af x , Af 2 and AP t ie-line for Step change in 
demand of 1% is applied at t=5s in hydro area shown in fig. 6(a)- 
(c) and fig.7(a)-(c) shows the responses of Af x , Af 2 and 

APtfe-ime f° r same Step load perturbation in thermal area. 

Examining these responses, it is clearly seen that the incremental 
changes in tie-line power is settling around 40s with PSO 
optimized I controller and around 90s for both I and PI controller 
when 1% SLP in hydro area. And better system performance in 
terms of minimum settling times in frequency and tie-line power 
deviations has achieved with PSO optimized I controller 
compared to conventional I and PI controllers. 

The peak overshoot and settling time are less in PSO 
optimized I controller and more in the conventional I controller 
for changes in the load in both hydro and thermal areas. 
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(a) 




(b) 




(c) 



Fig. 6 Comparison of integral, PI and I controller using PSO 
responses with 1% SLP in hydro area, (a) Af x , (b) Af 2 and (c) 



AP 



tie -lie 




Fig. 7 Comparison of integral, PI and I controller using PSO 
responses with 1% SEP in thermal area, (a) Afj , (b) Af 2 and (c) 



Table. 1 Comparison of time response specifications of frequency 
for single area AEFC 



S.n 

0 . 


Parameters 


Change in 
frequency 

(I) 


Change in 
frequency 
(PI) ' 


Change in 
frequency 
(PID) 


1 


Peak overshoot 


52% 


32.5 % 


9% 


2 


Settling time 


10 S 


18 S 


19 S 


3 


Delay time 


2.2 S 


2.2 S 


2.2 S 


4 


Rise time 


2.4 S 


2.3 S 


2.25 S 


5 


Peak time 


2.8 S 


2.5 S 


2.3 S 



Hence, it can be concluded that the optimized control 
strategy, PSO optimized I controller, provides a better control 
under wide variations in the loading condition of the system. 




(a) 



Conclusion 

In electricity power industry, there is an ongoing need for 
efficient and effective EFC techniques to counter the ever- 
increasing complexity of large-scale power systems and 
robustness against parameter uncertainties as well as plant/model 
mismatch and external load change. In this paper, a single area 
thermal power system is considered at the first instance with 
conventional controllers, out of these the PID controller provides 
superior results compared to I and PI controllers. And a widely 
used two-area hydro thermal system problems with I, PI and PSO 
optimized I controller are presented. Among these, with PSO has 
got the best dynamic response of frequency and tie-line powers. 
Future work may be carried out by including the study of the 
effect of changes in loading and system parameters with these 
controllers in the interconnected multi (Hydro-thermal-PV) area 
Automatic Foad Frequency controls. 
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APPENDIX 

Nominal parameters of Interconnected Hydro -thermal system 
investigated: 

f = 60Hz Rj = R 2 = 2.4Hz / p.u.MW 
T g - 0.08s = 200MW 

T r = 10s K r = 0.5 T t = 0.3s 

K d = 4.0 K p = 1.0 K; =5.0 

P rl = P r2 = 2000 MW T w = 1.0s 

Pi = P 2 = 8.33 * 10 ” 3 p.u.MW /Hz 
a 12 =-l T 12 =0.086 p.u.MW/ rad 

K pl = K p2 = 120 Hz / p.u.MW 
T p i = T p2 = 20s 
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